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Abstract 

A method is presented for calculating the Lie point symmetries of 
a scalar difference equation on a two-dimensional lattice. The sym- 
metry transformations act on the equations and on the lattice. They 
take solutions into solutions and can be used to perform symmetry 
reduction. The method generalizes one presented in a recent publica- 
tion for the case of ordinary difference equations. In turn, it can easily 
be generalized to difference systems involving an arbitrary number of 
dependent and independent variables. 



1 Introduction 

A recent article [1] was devoted to Lie point symmetries, acting on ordinary 
difference equations and lattices, while leaving their set of solutions invariant. 
The purpose of this article is to extend the previously obtained methods and 
results to the case of partial difference equations, i.e. equations involving 
more than one independent variable. 
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Algebraic techniques, making use of Lie groups and Lie algebras, have 
proved themselves to be extremely useful in the theory of differential equa- 
tions [2]. 

When applying similar algebraic methods to difference equations, several 
decisions have to be made. 

The first decision is a conceptual one. One can consider difference equa- 
tions and lattices as given objects to be studied. The aim then is to provide 
tools for solving these equations, simplifying the equations, classifying equa- 
tions and their solutions, and identifying integrable, or linearizable difference 
equations [1,3-26]. Alternatively, one can consider difference equations and 
the lattices on which they are defined, to be auxiliary objects. They are 
introduced in order to study solutions of differential equations, numerically 
or otherwise. The question to be asked in this is: how does one discretize a 
differential equation, while preserving its symmetry properties [27-31]. 

In this article we take the first point of view: the equation and the lattice 
are a priori given. The next decision to be made is a technical one: which 
aspect of symmetry to pursue. For differential equations one can look for 
point symmetries, or generalized ones. When restricting to point symmetries, 
and constructing the Lie algebra of the symmetry group, one can use vector 
fields acting on dependent and independent variables. Alternatively and 
equivalently, one can use evolutionary vector fields, acting only on dependent 
variables. For difference equations, these two approaches are in general not 
equivalent and may lead to different results, both of them correct and useful. 

Several aspects of symmetry for discrete equations were pursued in earlier 
articles by two of the present authors (D.L. and P.W.) and collaborators. The 
"intrinsic method" which provides, in an algorithmic way, all purely point 
symmetries of a given differential - difference equation on a given uniform 
fixed lattice was introduced in [4]. This was complemented by the "differen- 
tial equations method" in [5]. In addition to point symmetries the differential 
equation method provides a class of generalized symmetries. It was pointed 
out that in many cases the two methods provide the same result, i.e. all sym- 
metries are point ones. The two methods were successfully applied to many 
specific problems [5,7,12,13,16]. The advantage of these two approaches are 
their simplicity, their algorithmic character, and their close analogy to sym- 
metries of differential equations. Their disadvantage is that many interesting 
symmetries, like rotations among discrete variables, are lost in this approach. 

A complementary approach was first developed for linear difference equa- 
tions [8, 19], again given on fixed uniform lattices. It was formulated in terms 



2 



of linear difference operators, commuting with the linear operator defining 
the original difference equation. This approach provides a large number of 
symmetries and the symmetry algebras of the discrete equations and their 
continuous limits are actually isomorphic. The symmetries of the difference 
equations are not point ones: they act at many points of the lattice. They do 
however provide flows that commute with the flow determined by the original 
equation and can thus be used to obtain solutions. 

This aspect of commuting flows has been adapted to nonlinear difference 
and differential-difference equations [9-11,15,17]. The equations are defined 
on a fixed and uniform lattice. Generalized symmetries are considered to- 
gether with point ones and some of the generalized symmetries reduce to 
point ones in the continuous limit. The methods for finding these general- 
ized symmetries rely on either linearizability, as in the case of the discrete 
Burgers equation [9] , or on integrability (the existence of a Lax pair) as in the 
case of the Toda hierarchy [10, 11, 15], or the discrete nonlinear Schrodinger 
equation [17, 18]. 

This symmetry approach is powerful whenever it is applicable. Together 
with point and generalized symmetries it provides Backlund transformations 
as a composition of infinitely many higher symmetry transformations. This 
aspect has been explored in detail for the Toda lattice [15]. We emphasize 
that Backlund transformations for difference equations, just as for differential 
ones, are not obtained directly as Lie symmetries (not even as generalized 
ones). 

Each of the above methods has its own merits and will be further devel- 
oped in the future. 

In this article we take the same point of view as in our recent article [1]. 
We consider point symmetries only and use the formalism of vector fields 
acting on all variables, dependent and independent ones. In [1] we considered 
only one discretely changing variable. The lattice was not fixed. Instead it 
was given by a further difference equation. Point symmetries act on the 
entire difference system: the equation and the lattice. The lattice is not 
necessarily uniform and we explored the effect of choosing different types of 
lattices. The idea of using transforming lattices is due to Dorodnitsyn and 
coworkers [27-31]. We differ from them in one crucial aspect. They start 
from a given symmetry group and construct invariant difference schemes 
for a given group. We, on the other hand, start from a given difference 
scheme and find its Lie point symmetry group. Previously this was done 
for the case of one independent variable. In this article we generalize to 
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the multidimensional case. The generalization is by no means trivial. The 
lattice is given by N 2 equations, where N is the number of independent 
variables, all of them varying discretely. Transformations of continuously 
varying independent variables, if present, are also taken into account. 

We stress that the approach of this article complements those of previous 
ones. The results of [4] and [5] are obtained if we chose a special form of 
the lattice (e.g. x m+ i — x m = h in the case of one independent variable, 
where h is a fixed, nontransforming constant). We purpusely avoid any use 
of integrability. Like Lie theory for differential equations, this approach is 
applicable to arbitrary differential systems, integrable or not. 

A general formalism for determining the symmetry algebra is presented 
in Section 2. It generalizes the algorithm presented earlier [1] for ordinary 
difference equations to the case of several independent variables. In Section 3 
we apply the algorithm to a discrete linear heat equation which we consider 
on several different lattices, each providing its own symmetries. Section 4 is 
devoted to difference equations on lattices that are invariant under Lorentz 
transformations. In Section 5 we discuss two different discrete Burgers equa- 
tions, one linearizable, the other not. The lattices are the same in both cases, 
the symmetry algebras turn out to be different. Section 6 treats symmetries 
of differential- difference equations, i.e. equations involving both discrete and 
continuous variables. Some conclusions are drawn in the final Section 7. 

2 General symmetry formalism 
2.1 The difference scheme 

For clarity and brevity, let us consider one scalar equation for a continuous 
function of two (continuous) variables: u = u(x,t). A lattice will be a 
set of points Pj, lying in the plane M 2 and stretching in all directions with 
no boundaries. The points Pj in M 2 will be labeled by two discrete labels 
P m , n - The Cartesian coordinates of the point P m , n will be {x m , n ,t m>v ) with 
— oo < m < oo , — oo < n < oo (we are of course not obliged to use Cartesian 
coordinates). The value of the dependent variable in the point P„^ n will be 
denoted u m n '^(^x rn n ,t rn n J. 

A difference scheme will be a set of equations relating the values of 
{x,t, u} in a finite number of points. We start with one 'reference point' 
P m , n and define a finite number of points P m+ ^ n+ j in the neighborhood of 
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P m ,n- They must lie on two different curves, intersecting in P m ^ n . Thus, the 
difference scheme will have the form 

Eg, ^ {*^m+i,n+j > ^m+i,n+j; ^m+ijTH-.? } J 1 — ^ — ^ ^ ^ 

-ii < i < «2 - ji < j < 32 h, z"2, ji, J2 e Z-°. 

The situation is illustrated on Figure 1. It corresponds to a lattice deter- 
mined by 6 points. Our convention is that x increases as m grows, t increases 
as n grows (i.e. x m+1)n -x m>n = hi > , t m>n+1 -t m>n = h 2 > 0). The scheme 
on Figure 1 could be used e.g. to approximate a differential equation of third 
order in x, second in t. 




x 

Figure 1: Points on a lattice 



Of the above five equations in (JTJ, four determine the lattice, one the 
difference equation. If a continuous limit exists, it is a partial differential 
equation in two variables. The four equations determining the lattice will 
reduce to identities (like = 0). 

The system (JTJ) must satisfy certain independence criteria. Starting from 
the reference point P m ,n and a given number of neighboring points, it must 
be possible to calculate the values of {x, t, u} in all points. This requires a 
minimum of five equations: to be able to calculate the (x, t) in two directions 
and u in all points. For instance, to be able to move upward and to the right 
along the curves passing through P m , n (with either m, or n fixed) we impose 
a condition on the Jacobian 

_ d(Ei,E 2 , E 3: E 4 , E 5 ) 



+ 0. (2) 
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As an example of difference scheme, let us consider the simplest and 
most standard lattice, namely a uniformly spaced orthogonal lattice and a 
difference equation approximating the linear heat equation on this lattice. 
Equations in this case are: 

2<m+l,n •Km,n h\ tm+l,n ^m.n (3) 

■Km,n+1 •^m,n ^m,n+l tm,n ^2 (4) 

V>m,n+1 ^m,n t^m+l,n 2w mn -|- U m —\ n . . 

h = w {) 

where hi and h 2 are constants. 

The example is simple and the lattice and the lattice equations can be 
solved explicitly to give 

Xm,n = h 1 m + x t m<n = h 2 n + t . (6) 

The usual choice is xq = to = , hi — Y12 — 1 and then x is simply 
identified with m, t with n. We need the more complicated two index notation 
to describe arbitrary lattices and to formulate the symmetry algorithm (see 
below). 

The example suffices to bring out several points: 

1. Four equations are needed to describe the lattice. 

2. Four points are needed for equations of second order in x, first in t. Only 
three figure in the lattice equation, namely P m+ i in , P m , n and P m , n +i- To 
get the fourth point, P m _i jn , we shift m down by one unit in equations 

PSD. 

3. The independence condition (T5]) is needed to be able to solve for x m+ i }U , 

tm+l,ri) •Km,n+li ^m,n+l and Um^n+i. 



2.2 Symmetries of the difference scheme 

We are interested in point transformations of the type 

x = F\(x, t,u) t = G\(x, t,u) u = H\(x, t, u) (7) 

where A is a group parameter, such that when (x, t, u) satisfy the system ([1]) 
then (x, t, u) satisfy the same system. The transformation acts on the entire 
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space (x, t, u), at least locally, i.e. in some neighborhood of the reference point 
P m n , including all points P m+i)n+ j figuring in equation That means that 
the same functions F, G and H determine the transformation of all points. 
The transformations (jTJ) are generated by the vector field 

X = £(x, t, u)d x + t(x, t, u)d t + (j)(x, t, u)d u . (8) 

We wish to find the symmetry algebra of the system (pQ), that is the Lie 
algebra of the local symmetry group of local point transformations. To do 
this we must prolong the action of the vector field X from the reference 
point (x mtn , t m , n , «m,n) to all points figuring in the system (pQ). Since the 
transformations are given by the same functions F, G and H at all points, 
the prolongation of the vector field (jSJ) is obtained simply by evaluating the 
functions £, r and at the corresponding points. 

In order words, we can write 

-t-0(x m ri , t m n , u m n )d Ur 

where the summation is over all points figuring in the system (TjQ). The 
invariance requirement is formulated in terms of the prolonged vector field 
as 

piXE a \ Eb=0 l<a,b<5. (10) 

Just as in the case of ordinary difference equations, we can turn equation 
(flUj) into an algorithm for determining the symmetries, i.e. the coefficients 
in vector field (jSj). 

The procedure is as follows: 

1. Use the original equations (pQ) and the Jacobian condition §2§ to express 
five independent quantities in terms of the other ones, e.g. 



^1 •Em+i2,n ^2 tm+i2,n ^3 ■^m,n+j 
^4 tm,n J rj2 ^5 ^m+«2>™+i2 



as 



-h < i < %i - 1 - h < 3 < h ~ !■ 



(11) 



(12) 
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2. Write the five equations (TTUj) explicitly and replace the quantities v a us- 
ing equation (fT2"j) . We obtain five functional equations for the functions 
£, r and 0, evaluated at different point of the lattice. Once the func- 
tions v a are substituted into these equations, each value of x^k, and 
Ui^ is independent. Moreover, it can only figure via the corresponding 
£,i,ki T i,k and 4>i t k (with the same values of % and k), via the functions v a , 
or explicitly via the functions E a . 

3. Assume that the dependence of £, r and <fi on their variables is analytic. 
Convert the obtained functional equations into a system of differential 
equations by differentiating with respect to the variables a;^,^ and 
Ui^- This provides an overdetermined system of linear partial differen- 
tial equations which we must solve. 

4. The solutions of the differential equations must be substituted back 
into the functional ones and these in turn must be solved. 

The above algorithm provides us with the function £(x, t, u), r(x, t, u) and 
4>(x,t,u) figuring in equation ([H]). The finite transformations of the (local) 
Lie symmetry group are obtained in the usual manner, by integrating the 
vector field (JHJ): 

§ = £(x,t,u) f x =r(x,t,u) § = <f>(x,i,u) 

(13) 

x | A =o = x t \x=Q =t u \ x =o = u. 



3 Discrete heat equation 

The heat equation in one-dimension 

u t = u xx (14) 

is invariant under a six- dimensional Lie group, corresponding to translations 
in x and t, dilations, Galilei transformations, multiplication of u by a constant 
and expansions. It is also invariant under an infinite dimensional pseudo- 
group, corresponding to the linear superposition principle. 

Symmetries of the discrete heat equation have been studied, using differ- 
ent methods and imposing different restrictions on the symmetries [8,19,27, 
28]. 

Here we will use the discrete heat equation to illustrate the methods of 
Section 2 and to show the influence of the choice of the lattice. 
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3.1 Fixed rectangular lattice 

The discrete heat equation and a fixed lattice were given in equation (jSJ) and 
(j3J), (jl]), respectively. Applying the operator (jUJ) to the lattice, we obtain 

£(*£m+l,n; £m+l,n> ^m+l,n) £(-£m,n; tm,ni "m,ti) (^^) 

The values w m+ i jn , M m , n +i, w mjTl are not related by equation dSJ) (since it 
also contains w m _i,„). Hence if we differentiate equations ffT5l) . (TX6]) . e.g. with 
respect to u m , n , we find that £ is independent of u. We have £ TO +i, n = t m ,n 
so equation ffT5]) implies that £ does not depend on x. Similarly, equation 
f fT6|) implies that £ does not depend on £. Hence £ is constant. Similarly, 
we obtain that r(x,t,u) is also constant. Applying the prolongation prX to 
equation fl5J) we obtain the functional equation 

0m,n+l — 0m, n = T^T^ (0m+l,n — 20 mj „ + m _i jn ) (17) 
with e.g. 4> m ,n = ^\^m,ni tm,ni ^m,n) ■ 

In m>n+ i we replace w m , n +i, using equation ([S]). We then differentiate 
with respect to w m +i )n and again with respect to w m _i in . We obtain 

0m, n y 4(-^m,n; ^m,n)'^ , m,n B{x m n} t m ,n)- (^^) 

Substituting ffl8l) into equation (IT7|) . using (jSJ) again and setting the coef- 
ficient of u m +i, n , u m -i,n, u m , n and 1 equal to zero separately we find that 
A must be constant and B must be a solution of equation (jSJ). Thus, the 
symmetry algebra of the heat equation on the lattice ([3]), (j3J) is given by 

P = <9 t A = 9* W = «d„ § = S(x,t)d u (19) 

with 5 a solution of the equation itself. Thus, the only symmetries are those 
due to the fact that the equation is linear and autonomous. 

3.2 Lattices invariant under dilations 

There are at least two ways of making the discrete heat equation invariant 
under dilations. 
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A) Five point lattice 

We replace the system of equations Q, (j4]) and ([5]) by 

Xm+l,n 2x m rl -|- X m -.\ n 2'm.n+l X mn (^0) 

tm+l,n ^m.n ^rra,n+l 2t mn -|- t m ^ n ^\ (21) 

Um,n+1 U mtn ^m+l,7J 2lt m n -|- U m — \ >n 



22) 



Applying prX of equation (|S]) to (1211 and substituting for x m+ i >ni t m+ i >n , t m n+ i 
and x mjn+ i from the equations ff20|) . (J2TJ) we obtain 

(23) 

~ ^V^Vra— l,ri; ^m— 1,71) y>m— l,n) 

£ (•^'771,71) 2t mj n ^m,n—li ^m,n+l) £("^m,7i.) ^m,ni V"m,n)- (^4) 

Since it m ,n+i and M m>n are independent a differentiation of ff24|) with re- 
spects to say u m _i jn (contained on the left hand side via w m , n +i) implies that 
£ does not depend on it. Differentiating (12~41) with respect to t TOi „_i we find 
that £ cannot depend on t either. Putting £ = £(x) into equation (1231) and 
taking the second derivative with respect to x m _i >n and x m „, we obtain that 
£ is linear in x. Similarly, invariance of equation (!2~T1) restricts the form of 
r(x,t,u). Finally the lattice (1201) . (12T1) is invariant under the transformation 
generated by X with 

£ = ax + (3 T = ^t + 5. (25) 
Now let us apply prX to equation (|2"2"|) . We obtain 



- > m,n+l 0m, n 0m+l,n 20 m n -|- m _x ra . x ^-m.+l,n 2lt mi n -|- U m —\^ 



(2a- 7 )- 



^771,71+1 ^771,71 (*^77l+l,7l -^771, 7l) (•^'771+1,71 '^771,71' J 

(26) 

Taking the second derivative <9 Um+1 n d Um _ 1 n of equation (j2"6"l) after us- 
ing the equation fT22|) to eliminate w m>n+ i, we find mj „ = v4 mjn (a;, t)u mtn + 
B mtn (x,t). Substituting back into equation ff26l) we obtain A m>n = A = 
const., and see that B mjn (x,t) must satisfy the original difference system. 
Moreover, we obtain the restriction 7 = 2a. 
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Finally, on the lattice (I2"0"j) . (l2"Tj) the heat equation (I22p has a symmetry 
algebra generated by the operators (1T9"]) and the additional dilation operator 

D = xd x + 2td t . (27) 

We mention that the lattice equations ( 1201) . ( 12T1) can be solved to give 
x = am + b , t = cn + d. At first glance this seems to coincide with the 
lattice ([6]). The difference is that in equation ([6]) hi and /12 are fixed constants. 
Here a, b, c and d are integration constants that can be chosen arbitrarily In 
particular, they can be dilated. Hence the additional dilational symmetry. 



B) A four point lattice 

We only need four points to write the discrete heat equation, so it makes 
sense to write a four point lattice. Let us define the lattice by the equations 

tm+l,n tm,n tm,n+l tm,n c(x m +\^ n X m ^n} 0. (29) 

On this lattice the discrete heat equation (122]) simplifies to 

Um,n+1 ^m,n c(^7ra+l,n 2\i mn -|- 1t m _i n ). (30) 

Applying the same method as above, we find that invariance of the lattice 
implies £ = Ax + B, r = 2 At + C . Invariance of equation (1301) then implies 
= Du + S(x, t) where A, B, C and D are constants and S(x, t) solves the 
discrete heat equation. Thus, the discrete heat equation on the four point 
lattice (}2~gj) , (1251) is invariant under the same group as on the five point lattice 

(EQD, (ED- 



3.3 Exponential lattice 

Let us now consider a lattice that is neither equally spaced, nor orthogonal, 
given by the equations 

■Em+l,n 2x m r! , -|- X m ^ n —\ X m)J1 -)-i (l -\~ c) X m ^ n (31) 
^m,n+l tm,n ^ ^m+l,n ^m,n (32) 
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-60- \ \^ ^ 

Figure 2: Variables (x,t) as functions of m and n for the lattice equations 
f l3Tj) . fl32|) . The parameters and the integration constants are, respectively, 
c = a/2, h = 1 and a = 7r, /3 = 0, to = 0. 



with c ^ 0, — 1. These equations can be solved and explicitly the lattice is 

t = hn + t x = (1 + c) n (am + (3) (33) 

where to, ot and (3 are integration constants. Thus while t grows by constant 
increments, x grows with increments which vary exponentially with time (see 
Figure 2). Numerically this type of lattice may be useful if we can solve the 
equation asymptotically for large values of t and are interested in the small 
t behavior. 



The heat equation on lattice (13 ip . (1311) can be written as 

^"m.n+l ^m.n ^m+l,n 2lt mjn -|- U m _\ n (34) 

h {x m +l,n >£m,n) 

Applying the symmetry algorithm to the lattice equations (13"T|) . (f3"2|) we find 
that the symmetry algebra is restricted to 

X = [a x + b (1 + c) t/h ] d x + T d t + (j)(x, t, u) d u , (35) 

where a, b and r are arbitrary constants (whereas c and h are constants 
determining the lattice). Invariance of the equation (13"4|) implies a = in 
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(I3"5|) and restricts <p(x, t, u) to reflect linearity of the equation and nothing 
more. The resulting symmetry algebra has a basis consisting of 

P 1 = (l + cf h d x P = d t W = ud u S = S(x,t)d u (36) 

where S(x, t) satisfies the heat equation. We see that the system is no longer 
invariant under space translations, or rather, that these 'translations' become 
time dependent and thus simulate a transformation to a moving frame. 



3.4 Galilei invariant lattice 

Let us now consider the following difference scheme 

— 7 2 — = T2 e (37) 

•Em+l,n 2^C m)J1 -|- X m —i tn (39) 
(*£m+l,n •Em,n)T~2 (*^m,n+l •^m,n)^~l C (^^) 

where r l5 r 2 and ( are fixed constants. 

The lattice equations can be solved and we obtain 

tm,n = T 1 m + T 2 n + t X m , n = OT X 171 + ^— ^ ~j U + X (41) 

where a, to and xq are integration constants. The corresponding lattice is 
equally spaced and in general, nonorthogonal (see Figure 3). Indeed, the 
coordinate curves, corresponding to m — const and n = const, respectively, 
are 

x - x = a(t-t )- §-n 

(42) 



X — Xq 



(t-t ) + ±-m. 



These are two families of straight lines, orthogonal only in the special case 
(a 2 + 1)tit 2 = a(. If we choose 

ot x t 2 - C = (43) 

then the second family of coordinate lines in equation (|42p is parallel to the 
x axis. 
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Figure 3: Variables (x,t) as functions of m and n for the lattice equations 
(EH}, (139|) . (140]) . The parameters and the integration constants are, respec- 
tively, 7~i = 1, T2 = 2, £ = 2 and cr = 1, x = 0, to = 0. 



Invariance of equation (!38|) implies that in the vector field we have t(x, t, u) 
a = const. From the invariance of equation (l39l) we obtain £ = A(t) x + -B(t) 
with 

A(t m+ i tn ) = A{t mtn ) B{t m+ i :n ) — 2B(t mjn ) + _B(t m _i in ) = 0. (44) 

Finally, invariance of equation ( 1401) implies A(t) = and B(t) — (3t + 7 
where /3 and 7 are constants. Now let us apply the prolonged vector field to 
equation (1371) . We obtain = Ru + S(x, t) where S(x, t) satisfies the system 
(l3T|) ..... (l40l) . The symmetry algebra is given by 

P = d t P 1 = d x B = td x W = ud u S = S(x,t)d u . (45) 

Thus, the system is Galilei invariant with Galilei transformation generated 
by the operator B. 

Let us now consider the continuous limit of the system (l3Tt) ..... (l4T)]) . We 
use the solution fT4Tl) of the lattice equations (1381) . (I39I) . (T40l) and for simplicity 
restrict the constants by imposing equation (1431) . We have, from equation 
(SU), (BSD 

tm,n+l tm,n 7~2 •^m,n+l •Em,n 

(46) 
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The continuous limit is obtained by pushing T\ -C 1, t 2 « 1, C ^ 1 an d 
expanding both sides of equation (1571) into a Taylor series, keeping only the 
lowest order terms. The LHS of equation (l37j) gives 



= u t + 0(r 2 ) 

and the RHS is given by 

("ra — 1 ,n 2"U mjn -|- U m — l,n) 
("C") [^(•^r?i,n "I - ^^1) ^m,n "I - Tl) 2 i u(x mjn , tm,n) ~\~ ^(•^m,n <2"Tl; ^m,n ^l)] 

= + f w^t + ^2«tt + C(ti)- 
The continuous limit of the system (1371) . .... (1401) is 

2 1 

«t = M a + -u xt + — m« cr^O. (47) 

(J 

The symmetry algebra of this equation, for any value of a, is isomorphic to 
that of the heat equation. In addition to the pseudo-group of the superposi- 
tion principle, we have 

P = d t D — x d x + 2t d t - \u d u - cx d t 

K = tx d x + t 2 d t — \(t + \x 2 )u d u — c(x 2 d x + xtd t — \xu d u ) 

(48) 

P 1 =d x + cd t W = ud u 

B = td x — \xu d u — c(x d x — 2t d t ) — c 2 x d t c= 1/a. 

The fact that the commutation relations do not depend on c suggest that 
equation (1471) could be transformed into the heat equation. This is indeed 
the case and it suffices to put 

c\(2+c 2 )x + ct] 

u(x.t) = e Wf^P wia.S) , 
v > v (49) 

a = x + ct (3 = (1 + c 2 )(t - cx) 
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to obtain 

Wfi = w aa . (50) 

Notice that while the difference equation (1371) on the lattice fl38l) . (139]) . fHOl 
is Galilei invariant, this invariance is realized in a different manner, than for 
the continuous limit ( 1471) . To see this, compare the operator B of equation 
with that of equation 



4 Lorentz invariant equations 

The partial differential equation 

Uxy = f(u) (51) 

is invariant under the inhomogeneous Lorentz group, with its Lie algebra 
realized as 

X-i = d x X 2 = d y L — yd x - xd y (52) 

(for any function f(u)). In equation flBTl) x and y are 'light cone' coordinates. 
In the continuous case we can return to the usual space-time coordinates 
z = x + y, t = x — y, in which we have 

u zz - u tt = f(u) (53) 

instead of equation (loTj) and the Lorentz group is generated by 

P = d t P 1 = d z L = td z + zd t . (54) 

Let us now consider a discrete system, namely 

Um+l ,n+l ^m,n+l Um+l.n ^m,n r / \ /rr\ 

jz ; _ r v ? . — — - i = -f ( 55 ) 

V. J 'm+l,n J - , m,n}\ym,n+L ym,n) 

%in+l,n 2x m rl -|- X m -.\ n %m,n+l •^m,n (<-*6) 
ym,n+l 2?/ m rl -|- ym,n—l £/m+l,n Z/m,n 0. (57) 

Applying the operator pr X (with t replaced by y) of equation ([9} to equations 
(PD, (J52D we obtain 

f = Ac + C r] = By + D. (58) 
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Requesting the invariance of equation (I55p we find that must be linear 

<p = a(x,y)u + @(x,y). (59) 
The remaining determining equations yield a = «o = const, and 

df d 2 f 
(A + B)—L- + (a u m , n + (3(x, y)) —f- = 0. (60) 

Thus, for any function / = f(u) we obtain the symmetries (152]) . just as 
in the continuous case (they correspond to B — —A, cto — [3 — 0). As in 
the continuous case, the symmetry algebra can be larger for special choices 
of the function f(u). Let us analyze these cases. 

a) Nonlinear interaction 

We have /" 7^ 0, hence (3 = (3q = const. The function must then satisfy 

(A + B-a )f + (a u + P)f' = 0. (61) 

For a ^ we take 

f = u P p 0, 1 (62) 

(we have dropped some inessential constants). The system fl55|) . fl56l) . (JBTjl 
is, in this case, invariant under a four-dimensional group generated by the 
algebra (|52p . complemented by dilation 

2 

D = xd x + yd y + ud u . (63) 

1 — p 

For a = 0, (3 7^ we have 

/ = e n . (64) 
The algebra is again four- dimensional with the additional dilation 

D = xd x + yd y - 2d u . (65) 

b) Linear interaction f(u) = u 

The only elements of the Lie algebra additional to ( 1521) are 

D = ud u S(f3) = (3d u (66) 

where (3 satisfies the system ( |55l) . ( |56l) . ( IBTl) with /(m) = m. The presence of 
D and 5"(/?) is just a consequence of linearity. 
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c) Constant interaction f(u) = 1 

The additional elements of the Lie algebra are again a consequence of linear- 
ity, namely 

L = xd x + yd y + 2ud u S = [Sx(x) + S 2 (y)]d u (67) 

where Si(x) and S 2 (y) are arbitrary (because Si(x) + S 2 (y) is the general 
solution of equation fl55l) with f{u) = on the lattice fl56l) . fl57|) ). 

To find a discretization of equation floU]) . invariant under the group cor- 
responding to flo^|) is more difficult and we will not go into that here. 

As stressed in the Introduction, the methods of this article can be applied 
to any difference system, but they provide only point symmetries. We could 
treat the integrable discrete Liouville and Sine- Gordon equations of Faddeev 
[32], or Hirota [33], but would not otain the generalized symmetries that 
are of interest. The correct formalism to use for these equations is that of 
Ref. [11]. 



5 Discrete Burgers equation 

The continuous Burgers equation is written as 

u t = u xx + 2uu X} (68) 

or in potential form as 

v t = v xx + vl u = v x . (69) 

We shall determine the symmetry groups of two different discrete Burgers 
equations, both on the same lattice. The lattice is one of those used above 
for the heat equation, namely the four point lattice ([28]) . ff29l . Each of the 
four lattice equations involves at most three points. Hence, for any difference 
equation on this lattice, involving all four points, the symmetry algebra will 
be realized by vector fields of the form (JH]) with 

f = Ax + B r = 2At + D (70) 

where A, B and D are constants (see section 3.2B). 
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5.1 Nonintegrable discrete potential Burgers equation 

An absolutely straightforward discretization of equation fl69l) on the lattice 

(im am is 



"m,n+l Wm,n U m j r \^ n 2,U mjn -\- U m —\^ n ^ ( U m j r \^ n Um,n j (71) 

Applying the usual symmetry algorithm, we find a four-dimensional sym- 
metry algebra 

Pi = d x P = d t D = xd x + 2td t W = d u . (72) 

5.2 A linearizable discrete Burgers equation 

A different discrete Burgers equation was proposed recently [9]. It is lin- 
earizable by a discrete version of the Cole-Hopf transformation. Using the 
notation of this article, we write the linearizable equation as 

(1 ~\~ h x U m ^n)\^U m j r 2^n 2,Um-\-l,n -|- Ura.n h x Um+l,n{Um+2.ri ^m,?iy 
^771,71+1 "m,nTC ~ ■ ~ ;: ■ ~ '-. 

1 ~r Cll x |W m +l,n ^m,n i ^xU J m,nU J m, J r \,n\ 

(73) 

tm+l,n tm,n ^m/ra+l •^m,n 0. 

In equation (1731) c is a constant, but h x is a variable, subject to dilations. 
The determining equation is obtained in the usual manner. It involves the 
function m>n at all points figuring in equation f!73|) . and also the constant 
A of equation flTOj) . The equation is too long to be included here, but is 
straightforward to obtain. The variable that we choose to eliminate using 
equation (IT3"|) is u mtTl+ i. Differentiating twice with respect to u m+2 , n we obtain 



d 2 0m,n+l du mt7l+ i _ 

duL„,, du m +o n dul 



V,n+1 ua m+2,n ua m+2,n 



(74) 



We differentiate (1741) with respect to u m , n and then, separately, with respect 
to M m _i >n . We obtain two equations that are compatible for c(l + c) 2 h x {l + 
h x u mtn ) = 0. Otherwise they imply that (ft is linear in m: = a(x,t)u + 
fi(x,t). We have c ^ 0, h x ^ 0, but the case c = — 1 must be considered 
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separately. We first introduce the expression for into the determining 
equation and obtain, after a lengthy computation (using MAPLE): a = —A, 
(3 = 0. For c = — 1 we proceed differently, but got the same result. Finally, 
the Lie point symmetry algebra of the system (1751) . (128|) . (1291) has the basis 

P = d t Pi = d x D = xd x + 2td t -ud u . (75) 

This result should be compared with the symmetry algebra of equation 
(175|) on a fixed constant lattice, found earlier [9]. The symmetry algebra 
found there was five-dimensional. It was inherited from the heat equation, 
via the discrete Cole-Hopf transformation. It was realized in a 'discrete 
evolutionary formalism' by flows, commuting with the flow given by the 
Burgers equation. The symmetries found there were higher symmetries, and 
cannot be realized in terms of the vector fields of the form considered in this 
article. 

6 Symmetries of differential- difference equa- 
tions 

6.1 General comments 

Symmetries of differential-difference equations were discussed in our previous 
article [1]. Here we shall put them into the context of partial difference 
equations and consider a further example. As in the case of multiple discrete 
variables, we will consistently consider the action of vector fields at points in 
the space of independent and dependent variables. To do this we introduce a 
discrete independent variable n (or several such variables) and a continuous 
independent variable a (or a vector variable a). A point in the space of 
independent variables will be P n , a) its coordinates {x„ iQ ,, z„ )Q } where both x 
and z can be vectors. The form of the lattice is specified by some relations 
between x rejQ ,, z n ^ a and Un,a = u{x n ^ a ^ z nj0l ). 

We shall not present the general formalism here, but restrict to the case 
of one discretely varying variable z = z n , — oo < n < oo and either one 
continuous (time) variable (£), or two continuous variables (x,y). 

For instance, a uniform lattice that is time independent can be given by 
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the relations 

z n+l,a ~ 2 z n,a + z n-l,a = (76) 

z n,a z n, a' C^) 

tn+l,a ~ tn,a = 0. (78) 

where a' is a different value of the continuous variable a. 

Conditions (J77[) . (1781 are rather natural. They state that time is the same 
at each point of the lattice and that the lattice does not evolve in time. They 
are however not obligatory. Similarly, equation (!76|) is not obligatory. The 
solution of equations (!76j) . (1781 is of course trivial, namely 

z n = hn + zq t = t(a) (79) 

and we can identify t and a (t = a, h and z are constants). 

The prolongation of a vector field acting on a differential-difference scheme 
on the lattice (ITS"]) ..... (1781) will have the form 

(80) 



+ ... 



where the dots signify terms acting on time derivatives of u. Since u n Q ,, u nt0l > 
and u n+ i tCt are all independent, equations (1771) and (!78|) imply 

C = C(^n) r = r(t). (81) 

On any lattice satisfying equation (!77j) . (1781 we can simplify notation and 
write 

X = ((z) d z + r{t) d t + 0(z, t, u) d u . (82) 

Similarly for an equation with one discretely varying independent variable 
z and two continuous ones (x, y) one can impose 



< 2 ,z n,ai,ce2 z n—l,ai,a2 (83) 



z n,a' 1 ,oi2 z n,ai,a2 (84) 

z n,ai,a' 2 z n,ai,a2 
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2<n+l,ai,a2 *^n,ai,a2 ^ (85) 

Vn+l,ai,a2 Dn,a\,ci2 



Invariance of the conditions (154j) and (155]) then implies that the vector fields 
realizing the symmetry algebra have the form 

X = ((z) d z + £(x, y) d x + r)(x, y) d y + <p{z, x, y, u) d u . (86) 

We can again simplify notation identifying x = a,\, y = cxi and solving 
to give z n = hn + z (h and z constant). 



6.2 Examples 

We shall consider here just one example that brings out the role of the lat- 
tice equations very clearly The example is Toda field theory, or the two- 
dimensional Toda lattice [13,34,35]. It is given by the equation 

u n ,xy = e""- 1 """ - e Un - u " +1 (87) 

with u n = u(z n ,x,y). 

On the lattice (1551) ..... (1551) we start with equation (155|) and have 

i i 
prX = £(x,y)d x + ri(x,y)d y + ( n +k(z) d Zn+k + ^ <f) n+k d Un+k + d Un>xy 

k=-i fc=-i 

(88) 

where is calculated in the same way as for differential equations [2]. 
Applying ( 1551) to equations ( 1531) and ( 1571) we find 

£ = £(x) r] = 7](y) ( n = Az n + B <f> n = (3 n (x,y, z n ) (89) 

and we still have two equation to solve, namely 

Pn+l-Pn + tx + Vy = (90) 
Prw = 0. (91) 



On the lattice (1531) . .... (155|) z n+ i and independent. Hence we can 

differentiate (19"U|) with respect to z n+ i and find that P n+ i is independent of 
z n+ i and hence of n. We thus find a symmetry algebra generated by 

P x = d x P 2 = d y L = xd y -yd x S = d z D = zd z 

(92) 

U(k) = k(x) d u V(h) = h(y)d u 
22 



where k(x) and h(y) are arbitrary smooth functions. Notice that S and D 
act only on the lattice and U(k) and V(h) generate gauge transformations, 
acting only on the dependent variables. 

If we change the lattice to a fixed, nontransforming one, i.e. replace fl83l) 

by 

^n— 1 ,0:1,02 ^n,ai,Q2 h (93) 

h = const, the situation changes dramatically. We loose the dilation D of 
equation (1921) . however z n+ i and z n are now related by equation (|93|) . The 
solution of equation (I90l . (|9"T|) in this case is 

A, = ^(£c + Vy) + Kx) + %). (94) 

On this fixed lattice the Toda field equations are conformally invariant and 
the invariance algebra is spanned by 

X(f) = f{x) d x + lf'{x) d u Y{g) = g(y) d y + ^g'{y) d u 

(95) 

U(k) = k(x)d u V(h) = h(y)d u S = d z . 

We see that giving more freedom to the lattice (three points z n+ \, z n , z n -\ 
instead of two) may lead to a reduction of the symmetry group, rather than 
to an enhancement. For the Toda field theory the reduction is a drastic one: 
the two arbitrary functions f(x) and g(y) reduce to / = ax + b, g = —ay + d, 
respectively (and only the element D is added to the symmetry algebra). 



7 Conclusions and future outlook 

The main conclusion is that we have presented an algorithm for determining 
the Lie point symmetry group of a difference system, i.e. a difference equa- 
tion and the lattice it is defined on. The algorithm provides us with all Lie 
point symmetries of the system. In Ref. [1] we considered only one discretely 
varying independent variable. In this article we concentrated on the case 
of two such variables. The case of an arbitrary number of dependent and 
independent variables is completely analogous though it obviously involves 
more cumbersome notations and lengthier calculations. The problem of find- 
ing the symmetry group is reduced to solving linear functional equations. 
In turn, these are converted into an overdetermined system of linear partial 
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difference equations, just as in the case of differential equations. The fact 
that the determining equations are linear, even if the the studied equations 
are nonlinear, is due to the infinitesimal approach. 

The symmetry algorithm can be computerized, just as it has been for 
differential equations. 

In previous articles (other than Ref. [1]) we considered only one discretely 
varying variable and a fixed (nontransforming) lattice [4-18]. The coefficients 
in the vector fields, realizing the symmetry algebra, depended on variables 
evaluated at more than one point of the lattice, possibly infinitely many 
ones. Thus, one obtained generalized symmetries together with point ones. 
For integrable equations, including linear and linearizable ones, the symmetry 
structure can be quite rich [8-11, 15, 17, 18]. In the continuous limit some of 
the generalized symmetries reduce to point ones [11,17,18] and the structure 
of the symmetry algebra changes. 

A detailed comparison of various symmetry methods is postponed to a 
future article. Applications of Lie point symmetries, as well as generalized 
symmetries, to the solution of difference equations, will be given elsewhere. 
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